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Abstract 

We construct a class of exact commensurate and incommensurate standing wave 
(SW) solutions in a piecewise smooth analogue of the discrete non- linear Schrodinger 
(DNLS) model and present their linear stability analysis. In the case of the com- 
mensurate SW solutions the analysis reduces to the eigenvalue problem of a transfer 
matrix depending parametrically on the eigenfrequency. The spectrum of eigenfre- 
quencies and the corresponding eigenmodes can thereby be determined exactly. The 
spatial periodicity of a commensurate SW implies that the eigenmodes are of the 
Bloch form, characterised by an even number of Floquet multipliers. The spectrum is 
made up of bands that, in general, include a number of transition points correspond- 
ing to changes in the disposition of the Floquet multipliers. The latter characterise 
the different band segments. An alternative characterisation of the segments is in 
terms of the Krein signatures associated with the eigenfrequencies. When one or 
more parameters characterising the SW solution is made to vary, one occasionally 
encounters collisions between the band-edges or the intra-band transition points 
and, depending on the the Krein signatures of the colliding bands or segments, the 
spectrum may stretch out in the complex plane, leading to the onset of instability. 
We elucidate the correlation between the disposition of Floquet multipliers and the 
Krein signatures, presenting two specific examples where the SW possesses a definite 
window of stability, as distinct from the SWs obtained close to the anticontinuous 
and linear limits of the DNLS model. 
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1 Introduction 



The question of enumerating and classifying the possible excitations in a dis- 
crete non- linear lattice (e.g., a Fermi-Pasta-Ulam(FPU) chain) is too complex 
to be even addressed in its entirety. However, important classes of excita- 
tions have already been identified and studied exhaustively, and their physical 
implications are in the process of being worked out. One such class consists 
of localised solutions, namely the discrete breathers with periodic or quasi- 
periodic time-variation, while a related class includes the travelling breather 
solutions. A third class of solutions can be described as non-linear spatially 
periodic travelling waves analogous to the phonon modes in linear lattices. The 
existence of such travelling waves has been established in non-linear discrete 
Klein-Gordon (NDKG) lattices as well as in FPU chains (0; 0) close to the 
linear limit (i.e., for small amplitudes) as also in the vicinity of the so called 
anti-continuous limit (small inter-site coupling). In the latter case, the travel- 
ling waves were constructed as an infinite array of breathers (multi-breathers) 
with an appropriate phase introduced at each breather site. 

More recently, commensurate and incommensurate standing wave (SW) solu- 
tions have also been investigated in the literature, once again close to the two 
limits referred to above (jl;0) • Close to the linear limit, small amplitude SW 
solutions in the NDKG chain have been obtained from the DNLS approxima- 
tion by referring to orbits in a small neighbourhood of a fixed point in a related 
area-preserving mapping. In the weak-coupling situation, on the other hand, 
one obtains the SW solutions by continuation from the anti-continuous limit of 
a solution with an appropriate coding sequence. In either of these approaches 
one can obtain SW solutions both commensurate and incommensurate with 
the underlying lattice. Even chaotic spatial structures are not ruled out. 

The question of stability of these standing wave solutions is a more complex 
one. In the case of the travelling waves, it appears that they are afflicted with 
a type of modulational instability, as a result of which they break up into 
an array of localised pulses. The SWs in general also happen to be linearly 
unstable close to the small amplitude limit. For the commensurate SWs the 
eigenfrequencies form bands on the real axis, but these bands have a tendency 
to overlap and extend into the complex plane leading to oscillatory instabil- 
ities for arbitrarily small amplitudes. The spectra for incommensurate SWs, 
on the other hand, are generaly cantor-like and require a more sophisticated 
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analysis (|3|;|4j), but here again one finds that the SW's with even arbitrarily 
small amplitudes are prone to oscillatory instability. 

In the present paper we present an exact construction of commensurate and 
incommensurate SW's in a piece-wise smooth (PWS) system resembling the 
DNLS model and perform an exact linear stability analysis for the commensu- 
rate SW's. The linearised evolution equation for an SW in the model depends 
essentially on two parameters - the breather frequency u> and a nonlinearity 
parameter 7 and, in contrast to breathers obtained perturbatively from the 
linear limit, one can identify well-defined regions in the tu-7 parameter space 
where the SW's are linearly stable. In other words, away from the linear and 
anti-continuous limits there do exist stable commensurate SW solutions - the 
main result of this paper. In addition, the bands of eigenfrequencies arising 
from the linearised evolution equation which we calculate in arriving at the 
stability result present interesting features. While in general there exist gaps 
between the edges of successive bands and the width of the band gaps varies 
with the parameters u, 7, instability does not necessarily occur through a col- 
lision of the band-edges (i.e., shortening of a band-gap to zero width). As we 
shall see, there may exist, in addition to the band-edges, a number of transition 
points (see below) within the bands and collisions of these intra-band transi- 
tion points may equally well account for the onset of instabilty. The model 
under consideration being explicitly solvable, the detailed structure of a band 
can be worked out, including the location of the band-edges and transition 
points as also the eigenmodes associated with the different band segments. 
Additionally, the different band segments can be characterised in terms of the 
Krein signatures (jl; |fl 0) associated with the eigenfrequencies, and the cor- 
relation between the various possible dispositions of Floquet multipliers (see 
below) and the Krein signatures can be looked at. This gives us insight into 
infinite dimensional linear Hamiltonian systems with Bloch type eigenfunc- 
tions in the context of the theory of Hamiltonian Hopf bifurcation 0; Eol ) 
in infinite dimensional systems (see also (j?])). 

In section 2 below, we present the piecewise smooth DNLS model and briefly 
recall results from a couple of earlier papers (flit E2I 1 relating to the single-site 
and multi-site breather solutions. Section 3 presents the exact commensu- 
rate and incommensurate SW solutions, obtained by extension of the multi- 
site breathers where the breathers form an infinitely extended periodic or 
quasiperiodic array. Section 4 is devoted to the formulation of the stabil- 
ity problem with special reference to the commensurate SW solutions where 
the problem of calculating the spectrum of eigenfrequencies is reduced to the 
consideration of the eigenvalue problem of a 4 x 4 complex transfer matrix 
depending parametrically on the eigenfrequency p of the original problem as 
also on the breather parameters A (related to the breather frequency u>) , 7. We 
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indicate how the bands, together with their edges and transition points, are 
determined by the parametric bifurcation of the transfer matrix. This section 
also elucidates the correlation between the Floquet multipliers and the Krein 
signatures associated with the different band segments in the context of the 
linearised variational equations which constitute an infinite dimensional linear 
Hamiltonian system, and their role in characterising the onset of instability. 
In section 5, we specialise to a couple of particular cases involving SW's with 
specified spatial periodicity and phase characteristics. In both these cases the 
onset of instability is associated with a collision among intra-band transition 
points, and we obtain explicitly the window of stability in the A-7 plane. This 
section also contains concluding remarks. 



2 Breather solutions in a DNLS-like model 



In (HQ we considered a piecewise smooth DNLS-like model given by 

.d^, + +Vn-l) + f(\lpn\)lpn = 0, (la) 

at 

with 

/(a0 = 7(l--)e(*-l), (lb) 
x 

where stands for the Heavyside step function and 7 is a parameter char- 
acterising the strength of non-linearity. The non-linear term also involves a 
threshold parameter which we have scaled to unity. With reference to this 
threshold parameter we characterise lattice sites as being either high or low de- 



pending on whether \tp n \ > 1 or \ip n \ < 1 respectively. As indicated in (|llt ll^ 
the model is qualitatively similar to the cubic DNLS model (though the latter 
does not explicitly involve a threshold parameter) and possesses breather so- 
lutions of similar features. In particular, we presented explicit construction of 
single site and two-site breathers as also an explicit stability analysis of these 
breathers. 



The single site breather solution is characterised by two parameters : a spatial 
decay rate A (related to the breather frequency as uj = — A — \) and the 
nonlinearity parameter 7, and is given by 

i(j n = (2a) 



0n= J, ! A |n| . (2b) 

7 + A- j 

The parameter 7 is to be larger than a certain threshold value (depending on A) 
for this breather solution to exist. The linearisation of (la), (lb) around (2a), 
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(2b) in the rotating frame was then shown to possess a spectrum consisting of 
a single band with extended eigenmodes and an isolated point corresponding 
to a localised mode (in addition to a trivial phase mode with eigenfrequency 
zero). As the parameters A and 7 are made to vary, the isolated point in the 
spectrum gets shifted while the band remains fixed in the frequency scale and, 
as a certain stability border in the A-7 plane is crossed, the isolated point in the 
spectrum approaches zero (the 'double- zero' configuration) with a consequent 
onset of stationary instability. Interestingly, there exists a certain critical value 
of the spatial decay rate (A c ) above which the breather becomes intrinsically 
unstable. 

In (jl2l ) the model was also shown to possess 2-site monochromatic breather 
solutions of the form (2a) with high sites located at, say, n = and n = N, 
the amplitude at the high sites being 



Here 5 is the phase difference between the two high sites, which can be either 
or 7r, and once again 7 is to be larger than a certain threshold value depending 
on A for the breather to exist. An exact linear stability analysis was performed 
for the two-site breathers. The in-phase breather (5 = 0) was found to be 
intrinsically unstable while, for anti-phase breather (5 = tt), one obtains a 
band of extended modes together with a pair of localised symmetric modes 
and another pair of localised antisymmetric modes for parameters (A, 7) away 
from a stability border (for given N). As the parameters are made to vary, the 
isolated eigenfrequencies get shifted and, at the stability border the breather 
is destabilised through a Krein collision involving the symmetric eigenmodes. 

The model (la), (lb) admits of many other breather solutions including, pos- 
sibly, spatially random ones. We present below a class of exact commensurate 
and incommensurate standing wave (SW) solutions in the model and then per- 
form their linear stability analysis, with special reference to the commensurate 



3 Commensurate and Incommensurate SW solutions 



We consider a class of breather solutions for which the high sites are located in 
a periodic array at n = 0, ±iV, ±2N, ±3iV, .... However, successive high sites 
are allowed a phase difference 5 that may or may not be a rational multiple 



b 




(3) 



7 (l-A 2JV )-(i-A)(l-A w e i5 )' 



SWs. 
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of 2-7T (cf. the 2-site breather solution where 5 is either or n). In between the 
high sites there occurs a succession of low sites where one has 

^ + (^+1 + ^-1) =0. (4) 
Considering monochromatic breathers of the form (2a), (4) gives 

4> n = a\ n + /3A" ra , (5) 

where A satisfies 

A + j = -co, (6) 

and where we choose without loss of generality |A| < 1. Note that \u\ > 
2, i.e., the breather frequency lies outside the linear phonon band. In the 
following we restrict ourselves to A > {to < —2) for the sake of concreteness 
(corresponding results are obtained for A < 0). 



To be more specific, we look for monochromatic SW solutions of the form 
(2a) where high sites are located in a periodic array on the lattice at, say, 
n — 0, ±iV, ±2N, . . . (N > 2). Each high site is flanked on either side by a low 
site and the magnitude \v/j n \ at each high site is the same, say, b. However, 
we allow for a phase difference 5 between successive high sites. We note from 
(la), (lb) that for a low site (i.e., one with \ip n \ < 1), 

(fin + <fi n+l + n _i = 0, (7) 

i.e., <{> n is of the form (5), where A satisfies (6) and a, (3 are appropriate 
constants. For a high site with \tp n \ — b > 1, on the other hand, 

wbe"" + n+1 + n _i = - 7 6e*»(l - i), (8) 

where 6 n denotes the phase at the site under consideration. Taking together 
all these requirements, one arrives at the following SW solution: 

/ 1 _ pit \N p iS _ \ N \ 

^N +m = be™ [^^r>r + T^w xN - m ) . ( 9 ) 

where the index n for a site has been expressed in the form 

n = kN + m (0 < m < N, k = 0, ±1, ±2, . . .), (10) 

with k,m integers determined uniquely by n (for given N). Note that, while 
the phase advance (S) is uniform between successive high sites (n = kN, n = 
(k + l)N), it is not necessarily so for arbitrarily chosen pairs of successive sites 
n, n + 1. 
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It is also important to note that the standing wave solution (2a), (9) is valid 
only if 7, u satisfy appropriate constraints corresponding to the consistency 
conditions 

|0o|>l, W<1, (11) 

which imply, in general, that 7 is to be larger than a certain threshold value 
71(A) depending on iV and 8 (as already mentioned, similar conditions also 



apply for single-site and two-site breathers obtained in ([lit Il2f )). Expressions 



of 71(A) for the particular cases N = 2, 8 = n, and also for N = 3, 8 = n, 
will be found below. 



Evidently, (9) represents a spatially periodic, or commensurate standing wave 
solution if 8 is a rational multiple of 2n, say, 8 = (~)27r (r, q, co-prime integers, 
q > 0), in which case both the amplitude and phase are repeated at intervals 
of N = qN lattice sites. On the other hand, if is irrational, (9) represents 
an incommensurate standing wave since the phase of tp n never repeats itself, 
but comes arbitrarily close to any preassigned value 8 for some n or other. 



In section (5), we consider for the sake of concreteness two particular com- 
mensurate solutions with 8 = ir and with N = 2, 3 respectively. The solution 
for N = 2, 8 = tt (N = 4) reads (up to translation by a lattice site) 



0o = -02 = b, (12a) 

01 = 03 = 0, (12b) 

n+ jV = 0n, (12c) 

where 

6=-^-, (12d) 

7 + CO 

and where 7 is to satisfy the consistency condition 

7 > 1iW\n=2, 6=n = -u. (12e) 



Recall that in this paper we restrict to w < —2, while similar results are 
obtained for uj > 2. In a similar manner, one has, for N = 3, 8 = tt (periodicity 
N = 6) 



0o = -03 = b, 

01 = -02 = -04 = 05 = Z , 

1 — UJ 

4> n+ N = 0n, 



(13a) 
(13b) 
(13c) 
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where 



and where 



b= — 7(i - 1 a ; 2 ) — , ( i3d) 

7 (1 -A3) -1^(1 + A3)' 1 ; 



/ \\\ 1 + A + A 3 + A 4 
7> 7i(A)U=3,5=7r = • (13e) 

Having obtained the commensurate and incommensurate standing waves in 
our piecewise smooth DNLS-like model, we now turn to their linear stability 
analysis. 



4 Linear Stability Analysis 

4-1 The transfer matrix 

We consider a perturbation u n over the breather profile (p n of (2a), (9), (10) 
in the rotating frame and write (we drop the bar over <f> n ) 

^n=(0n + «n(*))e- iwt . (14) 

Using this in (la), (lb) and linearising, one obtains 



iu n + uu n + u n+l + u n _i = -7{w n - TTj-M-f 1 ~ if)} e ~ !)■ ( 15 ) 

This actually represents a linear Hamiltonian system, and we seek a harmonic 
solution of the form 

u n (t) = a n e ipt + b n e- tp * t , (16) 

where p (in general complex) is an eigenfrequency with eigenmodes specified 
by {a n }, {b n }. Substituting in (15), we obtain the equations satisfied by the 
a' n s, and b' n s : 

pa n = uia n + a n+ i + o n _i + -y{a n - T^rh(-r ~ %)} © (l&»| - 1), (17a) 



= - - bu - 7 R - - ^)} e (W - 1). (i7b) 
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These equations define a linear mapping 



a n 

a n-l 

b* 



"n+l 
b* 

u n+l 



Ar 



\K-iJ 

where, using (9), the matrix A n is obtained as 



p — uj — 7 (1 
1 



2b 



J_ p ~2ik5 
2b C 









(a \ 

a n-l 

K 
\K-i) 



7_ 2ikS 
'2ft C 











-P-w-7(1-s) - 1 
1 



(18) 



(19a) 



forn = A;A^ (k = 0, ±1, ±2, . . .), and 



A n = A 2 (p) 



V 



p-uj-l 
10 
-p-u -1 
10 



(19b) 



/ 



for n 7^ fciV. 



In other words, the mapping from high site to its next (low) site differs from 
the mapping from a low site to the next site, which may be a high or a low 
one. Since there are (N — 1) low sites in between two successive high sites, we 
obtain the transfer matrix from a high site at n — kN to the next high site 
n = (k + l)N as 

T k (p) = (A 2 (p)f~ 1 A[ k \p). (20) 
Thus, for instance, the mapping from the column vector 



X n = 



' ' On 


\ 




-1 


K 






-J 
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based at sites n = kN to the vector based at n = (k + r)N is 
X(k+r)N = T k+r ^i(p)....T k+1 (p) T k (p) X kN . 



(21) 



The criterion for the determination of the spectrum of eigenfrequencies can 
now be stated as follows : a complex number p is an eigenfrequency for the 
system (15) if there exists a vector X Q such that all the vectors X r ^ (r = 
0, ±1, ± 2, . . .) obtained from X by successive applications of T k (p) and 
T k ~ l (p) (with appropriate k) remain finite for r — > ± oo. Notice that for an 
incommensurate SW solution, for which 5 is not a rational multiple of 2n, the 
sequence of matrices T k (p) (k = 0, ±1, ± 2, . . .) never repeats itself, while 



for a commensurate solution with 5 
involves 25 rather than 5), 



: )27r one has (refer to (19a) which 



T k+s (p)=T k (p), 



(22) 



where s = q or | according as g is odd or even (s — 1 for 5 — 0). In the latter 
situation the problem reduces completely to that of diagonalising the 4x4 

transfer matrix 



R(p)=T 8 - 1 ( P )...T 1 (p)T (p). 



(23) 



In particular, for 5 = or 7r one has to diagonalise 



R(p) = (AMf- 1 A 1 (p) 



(24a) 



where 



Ai{p) 



P-w-7(l-s) "I 



26 

1 
_ _x 

2b 



26 



o N 



_ p _^_ 7 (i_ J_) _i 

1 y 



(24b) 
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and A 2 (p) is given by (19b). It is easy to check from (19b) that 



A 2 (p) 



N-l 



Hi N -Hi- N 
Mi-Mi -1 



„ 1 N - 1 - t , 1 - N + 1 w JV - 2 -m_- iV+2 

Mi-Mi -1 











V2 N -V2- N 



M2 iV - 1 -M2- iV + 1 



M2-M2" 



M2 iV - 1 -M2_- iV + 1 ^2 iV - 2 -M2_- iV + 2 
M2-M2" 1 M2-M2 -1 



/ 



(24c) 



where 



A*i,2 = 2 [±P ~ w ~ V (P =F ~ 4 1- 



(24d) 



Further, one observes that the eigenvalue of the transfer matrix R(p) occur in 
reciprocal pairs (A, 4) since positive and negative directions along the lattice 
are equivalent (it is easy to check this explicitly for (24a)). 



Referring to (16) one observes that an arbitrarily chosen value of p (in gen- 
eral complex) will be an eigenefrequency of the system (15) if at least one 
pair of the eigenvalues of the transfer matrix lies on the unit circle in the 
complex plane. One can thus determine the entire eigenvalue spectrum and, 
in addition, the eigenmodes. Evidently, the spectrum as also the eigenmodes 
depends parametrically on A, 7 (for given N, 5). If, for some particular choice 
of A, 7, the entire spectrum of eigenfrequencies lies on the real axis then the 
corresponding SW solution given by (9, 10) is linearly stable. 



One notes from (17a), (17b) that if p is an eigenfrequency then so are p*, —p, —p* 
(actually (15) is a linear Hamiltonian system; see section 4.2). Further, it fol- 
lows from above that, for the commensurate SW's, the eigenmodes are of 
the Bloch form. It is also apparent that the spectrum in this case consists of 
bands (as we see in sec. 5, the bands can be explicitly obtained in the present 
model). Starting from any A, 7 for which the entire band spectrum is con- 
fined to the real axis, one can vary one of the parameters (or an appropriate 
function of these) such that, as some critical parameter value is crossed the 
spectrum acquires complex eigenfrequencies. This corresponds to the onset of 
instability. 
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In the following, we apply these principles in sections 4.3, 5, to see how the 
bands, together with the band-edges and the transition points (see below) can 
be explicitly calculated and calculate the exact stability threshold for a couple 
of concrete cases (N = 2, 5 = n , and N = 3, 5 = n). However, before that we 
digress briefly to explain two important characteristics of the eigenfrequencies 
and the eigenmodes, namely the Krein signatures and Floquet multipliers. 

4-2 Linear Hamiltonian systems with spatial periodicity : Krein signatures 
and Floquet multipliers 



The evolution of the M n 's (eq. (15)) can be obtained from the Hamiltonian 

H = Yji( U + 7) Kl 2 + Un+lK + U nK+l + 
n 

+ 2hb HM "' 2 + 2|^ (0n2<2 + U «^* 2 V 9 (1^1 " (25) 
Note that, clS db result, the antisymmetric product 

2H = -i ^{u n u* n - u*u n }, (26) 
is conserved. Assuming a harmonic time dependence 

u n = a n e ipt + b n e- vpH , (27) 

we find 

-z£K< - u* n u n ) = (p + p*) {E(l«n| 2 - |^| 2 )} e l(p - pt)t - (28) 



Thus, in order that the left hand side be conserved we must have 

E(kl 2 -N 2 ) = o, (29) 

n 

if p 7^ p* . On the other hand, for a real p, the signature of the quantity 
XXI a n| 2 — l^n| 2 ) which is proportional to the energy evaluated for the solution 
(27), is termed the Krein Signature (see (Q)) for frequency p corresponding 
to the eigenmode specified by {a n }, {b n }. Denoting the eigenmode by the 
symbol x we express the Krein signature as K(p;x)- We have seen above 
that K(p; \) is zero if Im(p) 7^ 0. If the eigenfrequency p is degenerate with 
eigenmodes xi; • • -?Xs then for each eigenvector we have a Krein signature 
K(P'i Xt) (i = 1, ■ ■ ■ s). More generally, we can talk of the Krein signature for a 
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subspace, say, V, spanned by a given subset of the set of eigenvectors xu ■■■■Xs, 
and call it K(p;V). It is defined to be '+', '— ', or if the reduced energy 
i(p^p«) restricted to the subspace V is positive definite, negative definite, or 
indefinite (0). Finally, if V is taken to be the span of all the eigenvectors 
Xi, ■ ■ ■ , Xs, then K(p; V) is termed the Krein signature for the eigenfrequency 
p and denoted simply as K(p). 



The relevance of the Krein signature lies in the fact that it characterises crit- 
ical changes in the structure of the spectrum of eigenfrequencies for the given 
quadratic Hamiltonian H as a function of parameters through which H can 
be made to change continuously. In particular, it indicates the onset of bi- 
furcations wherein the the spectrum acquires eigenfrequencies with non-zero 
imaginary parts. For instance, we have the important result (0;0) that if v be 
a relevant parameter in H such that two real isolated eigenfrequencies P\{y), 
P2{v) with opposite Krein signatures collide at p as v approaches a value u, 
say, from below, then after the collision as 5 = (y — v) becomes positive the 
two eigenfrequencies spread out in the complex plane as p + ± ^2(8), 
where e\ and 62 are functions tending to zero as 5 — > ( the converse of this 
statement is also true, see (0)). Such collisions are responsible for the so called 
Hamiltonian Hopf bifurcation in non-linear Hamiltonian systems - a subject 
widely discussed in the literature (see (0 H; 0; E3) an d references therein). 



While the eigenfrequencies and the eigenmodes of a linear Hamiltonian system 
are characterised by their Krein signatures, much more can be said about 
the eigenmodes if the system happens to possess spatial periodicity. Thus, in 
(25), the coefficients of the linear system are spatially periodic because of the 
periodicity of the breather profile function (<f) n ) for the standing wave solution 
under consideration, and hence the eigenmodes have a Bloch-like structure, 
characterised by appropriate Bloch wave numbers. Thus, the coefficients a n 
(as also b n ) in (27) are of the form 

a n = r n A n , (30) 

where A n = A n+ fj (N is the spatial period ; N may be different from N, 
the spatial period of the commensurate SW; e.g., for 5 = it, N = 2N, but 
N = N.) , and r is the Bloch multiplier. One can simplify by looking at the 
values of a n at n — k N (k = 0, ±1, ±2, . . .) so that, writing a k ^ = dk we have 

~a k = AX k , (31) 

where A is a constant and A = T N is the Floquet multiplier characterising the 
eigenmode. As a consequence of the fact that both the two spatial directions 
in the lattice are equivalent, Floquet multipliers necessarily occur in pairs 
like (A, \). Choosing an arbitrarily specified frequency p (in general complex), 
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the solution for a n (also for b n ) can be written as a linear combination of 
special solutions of the form (30), (31) with Floquet multipliers, say, (Aj, j-) 
(i — 1,2,...). If at least one of these pairs lies on the unit circle (i.e., | Aj| = 1 
for some i) then the chosen value of p qualifies as an eigenfrequency, with these 
specified A/s determining the eigenmodes Xj (j — 1 ; 2, . . . , s, say) belonging to 
p. If there be r pairs of multipliers lying on the unit circle, then the maximum 
possible number of eigenmodes is s = 2r. However, s can be less than 2r 
depending on degeneracies among the X's. 

When a relevant parameter v characterising the Hamiltonian is made to vary, 
the eigenfrequency p, in general, moves continuously in the complex plane and, 
at the same time, the Floquet multipliers also move in the complex A- plane. 
In the process, certain critical values of v may arise when pairs of eigenvalues 
either leave or enter into the unit circle in the A-plane. The range of values 
covered by p for which at least one A-pair resides on the unit circle constitutes 
a band while inside the band, one may have transition points corresponding to 
an increase or decrease in the numbers of such pairs as mentioned above. The 
transition points divide a band into segments and at certain special values of 
the parameter u, a segment may shrink to zero width through a collision of a 
pair of transition points. As we shall see below, it is precisely at these special 
values that the Krein signature undergoes a change indicating, possibly, the 
onset of an instability. 

In summary, the Krein signatures associated with the eigenfrequencies form- 
ing a band in an infinite dimensional linear Hamiltonian system with Bloch 
type eigenmodes are intimately connected with the Floquet multipliers charac- 
terising these eigenmodes. In particular, collisions of transition points within 
bands (as also collisions between band-edges) are associated with changes in 
the Krein signature signifying a possible onset of instability. 



4-3 Determination of the band-edges and transition points. 

Referring to sections 4.1, and 4.2, we observe that the eigenvalues of the 
transfer matrix R for a commensurate SW solution, are precisely the Flo- 
quet multipliers for the problem. There are two pairs of these multipliers for 
any (complex) p, of which at least one pair has to be on the unit circle for 
p to qualify as an eigenfrequency. The multipliers, in general are of the form 
(Ai, A2, ^). However, for real p, the transfer matrix R(p) is real, and 
a multiplier A implies others of the form j-, A*, 4r. Assuming that, for any 
given real p all four multipliers lie on the unit circle (fig. 1, Al) one observes 
that, as the eigenfrequency is made to vary along the real line within a band, 
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the two pairs move on the unit circle till one pair of multipliers collide at +1 
or —1 (fig. 1, A2, A3), after which that pair moves off the unit circle along 
the real line, (fig. 1, Bl, B2) while another pair of multipliers continues to lie 
on the unit circle. The degenerate configurations A2, A3 of fig.l thus mark 
transition points in the band of eigenfrequencies. At these transition points 
one has, respectively, either 

P + 2 - 2T = 0, (32a) 

or 

P + 2 + 2T = 0, (32b) 
where P{p) and T(p) are respectively given by 

T(p) = J2K (33a) 
P{p)= y £\ i X j , (33b) 

i<j 

in terms of the eigenvalues Aj (i = 1,2,3,4) of R{p). Away from a transition 
point, a possible scenario consists of the pair of multipliers remaining on the 
unit circle moving along the latter till these, in turn, become degenerate at +1 
or —1 (fig. 1, CI, C2, C3, C4) and then cease to lie on the unit circle (fig.l., 
Dl, D2, D3). These correspond to band-edges i.e., points on the real ]?-axis 
where the band terminates. At these band-edges one again has either (32a) or 
(32b), but now lying in a different region in the P-T plane (we do not enter 
here into a characterisation of these regions). A band on the real p-axis can 
also terminate as both pairs of multipliers on the unit circle collide to form 
a degenerate Krein configuration (fig.l, El), and then all four leave the unit 
circle (fig. 1, E2). Such a band-edge is characterised by the relation 

4(P-2) = T 2 , |T|<4, (34) 

between T(p) and P{p) defined above. 

Having identified the conditions (in the following we refer to Eqns.(32a),(32b), 
(34) as conditions I, 11 and III respectively) characterising the band-edges and 
transition points, we can then completely determine the bands lying on the 
real p-axis by referring to the secular equation of R(p) given by (24a)-(24c). We 
carry out this programme in the following two subsections, where we find that 
a band may have quite a complex structure involving several transition points 
dividing the band into segments. Each segment is characterised by a definite 
disposition of multipliers given by Al, Bl, or B2 of fig. 1. For an eigenfrequency 
p within a segment with disposition Al, there are four independent eigenmodes 
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Fig. 1. Schematic diagram depicting a number of possible dispositions of Floquet 
multipliers with reference to the unit circle, for a real eigenfrequency p 

Xi, X2, X3, X4, each of the form 



a kN 

a (k-l)N 

h* - 
kN 

b* 



\ 



D ik<f> 



P 
7 



(35) 



\ (k-l)N j 

where a, /3, 7, 5 can be explicitly calculated. One thus has the Krein signatures 

K(p, Xi) = K(p, X2) = o"i (say), (36a) 



K(p, Xs) = K(p, Xi) = 02 (say), (36b) 
and the Krein signature K(p) associated with the four-dim space made up of 
Xii X2, X3-> Xi is '+') ' — \ or according as G\ , (T 2 are both positive, both 
negative or of opposite signs. For configurations Bl, B2 on the other hand, 
there are only two independent eigenvectors X11X2 for any given p and one 
has in this case, 

K(p, xi) = K(p, X 2) = K{p) = a, (say) (37) 
where a can be either '+' or ' — '. 



A transition point or a band-edge, however, corresponds to a degenerate config- 
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uration of the multipliers and presents a special case. Considering a transition 
point, for example, with eigenfrequency p (say) with multipliers as in A2 of 
fig. 1, the multipliers are e ± *^(say), 1, 1 . The transfer matrix R(p) possesses 
only one eigenvector for eigenvalue 1, while an independent generalised eigen- 
vector can also be obtained. Correspondingly, one obtains an eigenmode \ of 
the form 

/ 



(38) 



a kN 






a (k~l)N 




P 


b* - 

u kN 




7 


^(fe-l)JV J 







and the Krein signature of this mode is 



K(p, X ) = s l gn(\a\ 2 + \(3\ 2 -\ 1 \ 2 -\5\ 2 ). (39) 

There exists a second eigenmode x' associated with the same multiplier +1 
that can be constructed using the eigenvector and generalised eigenvector of 
R(p) referred to above, but it diverges linearly with site index n. The Krein 
signature for this eigenmode is, in general, not defined. It can, however, be de- 
fined for a p corresponding to a degenarate transition point (i.e., one where two 
transition points in a band coalesce). One finds that, for such a p, K(p, x) — 
and the Krein signature corresponding to the subspace spanned by Xi x' is 
also zero. 



As indicated in sec. 4.2, the coalescence of a pair of transition points within 
a band or of two band-edges may mark the onset of instability of the SW 
under consideration. The SW becomes unstable by way of the eigenspectrum 
acquiring a complex part. Thus, there exist complex p- values for which at 
least one pair of Floquet multipliers reside on the unit circle. One thereby has 
a band stretching out into the complex plane with non-zero imaginary part 
which terminates, once again, as any one of conditions (32a), (32b) is satisfied. 
As we have seen in sec 4.2, the Krein signature for a p lying within this part 
of the spectrum is zero. 



In the following section we present explicit results for SW's with N = 2, 5 = n 
and N = 3, 5 = n, exploring the detailed structure of the bands, and looking 
at the Krein signatures associated with the transition points that collide and 
lead to instabilities of the SW solutions. We thereby obtain the stability border 
in the A-7 parameter plane demarcating stable from unstable SW solutions. 
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5 Details of band structure 



5.1 N = 2,5 = Ti 



We now apply the above considerations to determine the detailed structure of 
the bands including the band-edges and transition points and to obtain the 
criterion for the onset of instability. We consider first the simplest of the SW 
solutions, namely those with N = 2, 5 = n. Here the high sites alternate with 
low sites, and successive high sites have a phase difference of ix (SW solutions 
with N = 2, 6 = are intrinsically unstable). 



As we have seen in section 3, such an SW solution is given by (12a)-(12e) 
(recall that we confine ourselves in this paper to A > 0, oj < —2). One then 
has the transfer matrix 

R(p)=A 2 (p)A 1 (p), (40) 

where A ± (p), A 2 (p) are given by (24b), (19b) respectively, and it is easy to 
set up the secular equation for R(p), obtaining 

Tip) =J2 X i = 2 P 2 + (^ + w 7 - 4), (41a) 
P{p) = J2 = P A ~ (4 + co 2 ) p 2 - 2(u 2 + u; 7 - 3) (41b) 



One can now locate the band-edges and transition points from conditions I, 
II and III (eqns. (32a), (32b), (34) respectively). Using the notation 

z = p 2 , (42) 

One observes that T and P are polynomials in z of degree 1 and 2 respectively, 
and finds, 

(Condition I) z = z lj2 = -{8 + uj 2 ± y / (a; 4 + 32w 2 + I6W7)}, (43a) 

2 



(Condition II) z — 0, u 2 (= z 3 , say), (43b) 

(Condition III) z = z^ = — ^| ^~2~^~) ' (43c) 
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We mention that the solution (43c) satisfies |T| < 4 (see (34)) only if 

7 < 7 < -3a;, (44a) 

where 

7 = i[-2(o; 2 -4)- yV + 64)]. (44b) 



Recall that the standing wave solution under consideration exists only for 
7 > —uj. In order to describe the stability characteristics of the SW solution 
we look at (43a)-(43c) for various fixed values of uj (uj < —2), and, for each uj, 
consider possible values of 7 > —uj. For any given uj, the solutions (43b) are 
non-negative irrespective of 7 and hence correspond to real eigenfrequencies. 
The existence of the eigenfrequency p = for all relevant uj, 7 corresponds 
to the fact that the SW under consideration is a time-periodic solution of an 
autonomous set of differential equations, and relates to a phase shift of the 
periodic solution. On the other hand, the solutions (43a) are both real and 
positive for 

uj 3 

-uj < 7 < -2uj - — . (45) 

3 

At 7 = 7 C = — 2uj — ^ , the two solutions coalesce and, for 7 > 7,, they become 
complex, corresponding to complex solutions for p. Finally, one finds that in 
the range — uj < 7 < j c , the solution (43c) subject to (44a), if it exists, is 
positive, leading once again to real eigenfrequencies p. 



These observation lead us to conclude that all the transition points and band- 
edges, and hence the entire eigenfrequency spectrum, lie on the real p-axis 
for any given uj (< —2) and 7 satisfying (45). Indeed, for such a (a;, 7) the 
spectrum in the right half plane consists of one or two bands with transition 
points in their interior. For the sake of illustration we assume that uj satisfies 
uj 2 < 8 and, for such a given uj, look at different values of 7 in the range (44a). 
In this case, z<± is positive and satisfies (34), and one has 

z 3 < z(= p 2 ) < z 4j (46a) 

where 




(46b) 



Considering, first, values of 7 in the range 

uj 3 

7<7< -2^-— (< -3cu), (47) 

the band of eigenfrequencies extends on the real line from p = to p = y/z~i 
(fig. 2(a); for the sake of convenience, we consider only the right half complex 
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Im p 




A 



Fig. 2. Schematic diagram of the spectrum for the standing wave solution with 
N = 2, S = it, for an appropriately chosen u and for three values of 7, respectively 
less than, equal to, and greater than the critical value, see text; dispositions of the 
Floquet multipliers for the different band segments are indicated in accordance with 
notation in fig. 1 



plane of eigenfrequencies; the part of the spectrum contained in the left half 
plane is obtained by reflection about the imaginary axis). 



The transition points in the band are at B (p — oj (= \fzs) ), and at C, C 
(Pi,2 = y/ Zi, 2), the latter two being the values of p obtained from condition 
I. As 7 is increased towards the value 7 C (for the chosen value of u satisfying 
u 2 < 8, the condition for Z3 to satisfy z 3 < z) the transition points C, C 
gradually come closer till, at 7 = 7 C they coalesce at A (p = p = \J\{ uj2 + 8)) 
(fig. 2(b)). As 7 is made to cross 7 C from below, the spectrum acquires a 
complex part and looks as in fig. 2(c). 



As mentioned above, Fig. 2(a,b,c), depicts schematically the detailed structure 
of the band for three different values of 7, namely 7 = 7 C — e, 7 = j c , 7 = 
7 C + e respectively, where e is sufficiently small. The dispositions of the Floquet 
multipliers in the different segments of the band are indicated in each case, 
in accordance with the notation of fig. 1 and fig. 3 (see below). One observes 
that at the critical value 7 = 7 C (fig. 2(b)) the segment CC of fig. 2(a) 
shrinks to the point A, i.e., the two transition points C,C corresponding to 
the solutions of condition I (configuration A2 of fig. 1) collide and, for 7 > 7 C 
(fig. 2(c)), these move off the real line into the complex plane. The spectrum 
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F2 



Fig. 3. Schematic diagram of a pair of possible dispositions of Floquet multipliers 
for a complex p. 

thereby includes the segment CAC stretching into the complex plane, and 
the standing wave solution becomes unstable. At the point A (fig. 2(c)), the 
disposition of multipliers is like Al of fig. 1, i.e, both pairs reside on the unit 
circle. For any p in the interior of the segments AC or AC, the disposition of 
multipliers is like that marked Fl in fig. 3 where one pair continues to reside 
on the unit circle while the other pair moves off into the complex A-plane (the 
stretches AC and AC in fig. 2(c) are not actually straight, but curved lines). 
At the edges C, C (fig. 2(c)), condition I is again satisfied (but with T, P 
lying in a region of the P-T plane different from the one corresponding to C 
or C in fig. 2(a)), and the disposition of multipliers is like F2 in Fig. 3. 



Let us consider, for 7 = 7 C — e, two points in the spectrum, one with p = pi — 8 
(slightly left of C) and the other with p = p 2 + 5 (slightly right of C) where 5 
is small. For each of these points there are two pairs of multipliers on the unit 
circle, of which one pair is close to +1. For p = p\— S, we denote the eigenmodes 
of the pair located away from +1 as xu X2 (they are complex conjugate to 
one another), and the remaining two eigenmodes ( for multipliers close to +1) 
as X3i X4 (again mutually complex conjugate). Similarly, for p = p 2 + 6 we 
denote the corresponding eigenmodes as x'n X21 X3 an d Xa- AH the eigenmodes 
have the general form 



02fe-l 

b*2k 
\b*2k-i J 



P 
7 



(48) 



for which the Krein signature is 

a = sign(\a\ 2 + \{3\ 2 



1 71 



l*| 2 )- 



(49) 



In the present case all the eigenmodes can be explicitly obtained and one finds 
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K(p 2 + S; x[) = K(p 2 + 5; X ' 2 ) = K( Pl - 5; X i) = K( Pl - 5; Xa ) = 



- (50a) 



K(p 2 +8;x' 3 ) = K(p 2 +5;x > d = -K(pi-S; X 3) = -K(p 1 -5;x4)= -. (50b) 



Thus, the Krein signatures associated with the eigenfrequencies pi — 5 and 
Pi + S are respectively 



These are also the signatures for all p in the interior of the segments C'B' and 
BC respectively. As regards the signatures for the transition points C and C 
(5 = 0) the situation is slightly different. In case of C (p = pi), for instance, 
there are again two eigenmodes Xi> X2 (complex conjugate of each other) of 
the form (48). For the multiplier +1, one eigenmode (X3, say) is again of the 
form (48) with (0 = 0), while there exits another eigenmode %4 that diverges 
linearly with site index n. The signature for xa is n °t wen defined while that 
for Xi-i X21 X3 are as in (50a,50b). Similar statements apply for the transition 
point C (p = p 2 ). 

At the collision (7 = j c , fig. 2(b)) where C and C coalesce at A, the Krein 
signatures for any eigenfrequency in the interior of AB' and BA are again as 
in (51a,51b) respectively. As regards the point A (p — p, 7 = j c ) the signature 
for the two eigenmodes Xi> Xi corresponding to the pair of multipliers located 
away from +1 is again '— ', while that for %3 is a = 0. The other eigenmode 
Xi for the multiplier +1 again diverges linearly; however, now its signature is 
well defined and one finds that the reduced energy restricted to the subspace 
spanned by X3, X4 is an indefinite quadratic form, so that one has, for 7 = 7 C , 



Summarising, we observe that for each of the colliding transition points C, 
C there exits a relevant Krein signature corresponding to ^3, x'3 respectively 
which are opposite and, at the collision, this relevant signature reduces to 
zero. Another way of characterising the collision is by referring to the Krein 
signatures for eigenfrequencies belonging to the interiors of the segments BC 
and B'C respectively. As seen from (50a, 50b), these are respectively and 
'— '. In other words, an open set of eigenfrequencies with signature a = is 
embedded within the spectrum both before and at the collision (see 



K( Pl -6) = 0, 
K(p 2 + 6) = -. 



(51a) 
(51b) 



K{p) = 0. 



(51c) 
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Fig. 4. Numerical integration of (la), (lb) to exhibit stability of N = 2, 5 = ir, 
SW solution (7 = 5.9765, slightly below the critical value) starting from the profile 
(as function of lattice site index n) at (a) t = up to r = 5000T (the new profile 
coincides with (a) ) ; (b) corresponds to t + |, and r + ^ , (c) tor + y, and r + ^ , 
and (d) tor + f . 

On the other side of the collision (7 > 7 C , fig. 2(c)), the Krein signatures for 
eigenfrequencies in the interiors of BA, AB' remain the same (i.e., and '— ' 
respectively), while that for the entire complex segment CAC one has a = 0. 



We thus have a complete description, both in terms of Floquet multipliers and 
Krein signatures, of the band structure and of the collision causing the onset 
of instability. As mentioned earlier, for every fixed A (and hence u = —A — 
there exits a critical value of 7 namely 7 = 7 c (w) = — 2cu — such that, as 7 
crosses 7 C from below, the SW breather becomes unstable. We show in fig. 4, 
fig. 5 this loss in stabilty for A = 0.5 (u = —2.5, 7 C ~ 5.9766) where we display 
the SW profile (eq.(2a), (12a)-(12e)), at t = as also the profile obtained by 
numerical integration of (la), (lb) up to a time r (say, see captions) for two 
values of 7, one just below and the other just above j c . One observes in these 
figures that the profile remains remarkably stable in fig. 4 for r as large as 
5000T, while in fig. 5, the profile quickly breaks up even at r = 1500T (T = 
breather time period). 

For each fixed value of A (0 < A < 1; uj < —2) the SW remains stable for 
a window of 7- values given by 71(A) < 7 < 7 C (A). As A is made to increase 
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Fig. 5. Numerical integration as in fig. 4, but with 7 = 5.9767, slightly above the 
critical value, where the profile obtained for r = 1500T is shown; it differs markedly 
from the breather profile, similar to (a) in fig. 4, at t = 0. 

from to 1 the width of the window decreases, i.e., for decreasing breather 
frequency, the onset of instability occurs at progressively lower values of the 
strength of nonlinearity 7 (fig. 6). 



It is interesting to compare, for the sake of completeness, a different type of 
collision of eigenfrequencies, namely one where the collision does not result in 
the appearance of complex frequencies. Consider, for instance, ui = ujq = — 3 
and 7 = 70 = y- With these values of oj and 7 (43a), (43b) give z\ = z 3 = 9. 
i.e., p = ±3. Fig. 7(a,b,c) shows schematically the eigenfrequency spectrum 
for uj = ujq = —3, and for 7 = 70 — e, 7 = 70, 7 = 7o + e respectively, where e 
is sufficiently small. 

The spectrum now consists of two bands, and the various parts of the spectrum 
marked by transition points and band-edges have once again been marked 
with labels indicating the respective dispositions of Floquet multipliers in 



with only one pair of multipliers residing on the unit circle for each p within 
this band (disposition B2 in fig.l) and the Krein signature for such a p is 
found to be '+'. We have thus labelled the band l B2+\ The spectrum for 
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Fig. 6. Window of stability of the N = 2,5 = it breather; the curve with the full line 
corresponds to the threshold value (71(A)) of 7 below which the SW solution does 
not exist, while the dotted curve corresponds to the stability border; for any fixed 
A, the vertical distance between the two curves depicts the 7-window of stability of 
the SW solution. 

co = —3, 7 = y — e contains a second band extending from B (p — —u = y/zi) 
to D(p = yfz±) through C (p — y/zT), the latter being a transition point 
satisfying (32a). For small e, the width of the segment BC is small and it is 
labelled 'B1-' to denote that the disposition of the Floquet multipliers for any 
p in the interior of the segment is Bl (fig. 1), and the Krein signature is ' — '. 
In a similar manner, the segment CD receives the label 'Al-'. 



As e — > i.e., 7 approaches the value y from below (note that this value is 
less than 7 C = 7.6875, the critical value corresponding to cu = — 3 at which one 
has the collision z\ = z 2 ) there occurs a collision between B and C, but now 
the two signatures involved are the same and there does not appear complex 
eigenfrequencies as a result of the collision, as seen in fig. 7(b,c). We observe 
from fig. 7(c) that the points B, C just pass through each other, and now the 
eigenfrequencies ^pz~2 form the band-edges. 



Observe that, for u> = —3, the critical collision, leading to instability occurs 
for 7 = 7c = 7.6875, when in fig. 7(c), the band-edges A and C approach each 
other. As seen from this figure, the Krein signatures involved are now opposite, 
and the spectrum acquires a complex part as 7 crosses 7 C from below. 



In summary, we have obtained the detailed band structure for iV = 2, 8 — ir, 
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Fig. 7. Schematic diagram of spectrum for N = 2, 5 = ir, SW solution similar to fig. 
2, but now with u = —3, and three values of 7, respectively less than, equal to, and 
greater than 

and have characterised the various bands and band segments in terms of the 
Floquet multipliers and Krein signatures, and have distinguished between two 
types of collisions in terms of the Krein signatures of the colliding bands or 
segments. 



5.2 Band Structure: N=3, 5 = n 



As another example of a collision of intraband transition points, we consider 
commensurate SW solutions with N = 3, 5 = 71 (solutions with = 3, 5 = 
are again intrinsically unstable) given by (13a)-(13e), where explicit formulae 
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Fig. 8. Schematic diagram, similar to fig. 2 and fig. 7, of the spectrum for the 
N = 3, 5 = it SW solution, with u = —2.5, and three values of 7, resp. smaller than, 
equal to, and larger than the critical value (see text). 



can once again be worked out. Using the notation (42) , we now have: 

2 



T 



-5a; + 



1-LO 



7)2; _(l_^(^ + 2)-(u; 2 -l) 7 , 



(52a) 



o , 9 2(1 + 7), 9 , a 4a; 3 o; 2 + a; + 2 N 

P = - z s + 2a; 2 + 8 - ^ /; )^ + + ; 9 + 2 7 — —)z + 

1 — a; 1 — a; 1 — a; 

(2(a; 3 -3a; + l) + 27(a; 2 -l)). (52b) 

We now have, in general, three eigenfrequencies (we once again restrict to 
Re(p) > 0) satisfying condition I (eq. (32a)), and three others satisfying con- 
dition II (eq.(32b)), of which one is the trivial solution p = 0, while the number 
of eigenfrequencies satisfying condition III depends on 7, to. We choose, for the 
sake of concreteness, A = 0.5 (a; = —2.5). Then two of the solutions to (32a) 
undergo a collision for 7 = 7 C = 4.75585 at p = 1.5750. We show schematically 
in fig. 8(a,b,c) the detailed band structures for three values of 7, namely 7 C — e, 
7c 7c + e respectively, where (e > 0) is small. 

For each of these three values of 7, we have two bands separated by a band 
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Fig. 9. Numerical integration of SW solution for N = 3, 5 = ir, with A = .5, 
7 = 4.75580 (slightly less than 7 C ); for explanation, see caption to fig. 4. 

gap, together with a number of transition points within one of the bands. 
As in fig. 7 we mark each of the band segments with a label indicating the 
disposition of Floquet multipliers as also the Krein signature characterising 
that segment. In fig. 8(a,b,c) the points B, C, E correspond to condition I, 
while O, A, F correspond to condition II, and D corresponds to condition III. 
The stretch from D to E is a band gap. One observes that, at 7 = 7 C , the 
transition points B and C within the first band collide and then stretch out 
into the complex plane (the segment from B to C in the complex plane in fig. 
8(c) is not actually straight but is slightly curved). The label in the segment 
from A to D in fig. 8(c) is meant to imply that the Krein signature for any 
p within this segment is (i.e., the reduced energy is an indefinite quadratic 
form on the space spanned by the eigenmodes) while the disposition of the 
Floquet multipliers is similar to Al in fig. 1. We note from fig. 8(b) that the 
two colliding segments bear opposite Krein signatures, as a result of which 
they stretch out into the complex plane on the other side of the collision. 

Figures 9, 10 show the results of numerical integration of an initial standing 
wave profile with A = 0.5, and with 7 =4.75580(< j c ) and 7=4.75590(> 7 C ) 
respectively, where we find that the SW solution indeed gets destabilised as 7 
crosses 7 C from below. 

Fig. 11 depicts the region of stability of the standing wave solution under 
consideration in the 7 — A parameter plane, where we once again find that, for 
decreasing breather frequency, the SW gets destabilised at progressively lower 
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Fig. 10. Numerical integration of SW solution for N = 3, <5 = ir, with A = .5, 
7 = 4.75590 (slightly greater than 7 C ); for explanation, see caption to fig. 5. 
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Fig. 11. Stability window for the N = 3, S = w, SW solution; for explanation see 
caption to fig. 6 

values of the non-linearity parameter 7. 



In conclusion, we observe that, away from the linear and the anti-continuous 
limits, there exist dynamically stable commensurate standing wave solutions 
in our DNLS-like model, in contrast to DNLS breathers close to the above 
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two limits. The frequency spectrum and eigenmodes of the linearised systems 
around these SW solutions present interesting features that can be described 
in terms of the dispositions of Floquet multipliers as also the Krein signatures 
associated with the various segments of the band spectrum. The onset of 
instability occurs not necessarily through collisions of band-edges, but also 
through collisions of intra-band transition points. 

In a future communication, we shall study exact commensurate and incom- 
mensurate SW solutions in a piecewise linear NDKG model and look into 
their stability problem, with special reference to incommensurate SW solu- 
tions. The possibility of multibreaher solutions involving an infinite array of 
breathers with random spatial distribution will also be investigated in these 
exactly solvable models. 
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